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ABSTRACT 

The question of possible small-scale dynamo action in the surface layers of the Sun is revisited with realistic 3-D MHD 
simulations. As in other MHD problems, dynamo action is found to be a sensitive function of the magnetic Prandtl 
number Pm = i^/ri; it disappears below a critical value Pc which is a function of the numerical resolution. At a grid 
spacing of 3.5 km. Pc based on the hyperdiffusivities implemented in the code (STAGGER) is « 1, increasing with 
increasing grid spacing. As in other settings, it remains uncertain whether small scale dynamo action is present in the 
astrophysical limit where Pm << 1 and magnetic Reynolds number Rm S> 1. The question is discussed in the context 
of the strong effect that external stray fields are observed to have in generating and maintaining dynamo action in 
other numerical and laboratory systems, and in connection with the type-II hypertransient behavior of dynamo action 
observed in the absence of such external fields. 
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1. Introduction 

In laboratory experiments and numerical simulations of hy¬ 
drodynamic turbulence, it is found that at sufficiently high 
Reynolds numbers the statistical properties of the flow be¬ 
come independent of the viscosity value (v). This obser¬ 
vation has been enshrined in the standard cascade picture 
for three-dimensional (but not two-dimensional) turbulence 
(Kolmogorov 1941). The addition of a magnetic field intro¬ 
duces a second dissipative mechanism, the magnetic diffu- 
sivity t]. By analogy with the hydrodynamics case, it may 
be assumed that a turbulent magnetic flow would also be in¬ 
sensitive to the value of rj, such that it would only affect the 
small scales in the magnetic field. Experience with numeri¬ 
cal magnetohydrodynamical (MHD) simulations has shown 
this assumption to be unexpectedly problematic, however. 
The behavior of turbulent MHD flows, including their bulk 
transport efficiency and the presence or absence of small- 
scale dynamo action, appears to vary with details of the 
problem studied, with the numerical methods used, and 
in particular appears to depend on the magnetic Prandtl 
number of the fluid. Pm = v jri = Rm/Re- 

The magnetic field at the surface of the Sun is domi¬ 
nated by the sunspot cycle, which is believed to have its 
source near the base of the convective envelope. At suffi¬ 
cient spatial resolution, a weak field of mixed polarity is 
also observed. Although it is not clear from the observa¬ 
tions whether this is a truly separate component and not 
some sort of waste product of the spot cycle, it suggests 
that a local small-scale turbulent dynamo process might be 
operating in the surface layers (Durney et al. 1993; Petrovay 
& Szakaly 1993). This suggestion has been addressed with 
a number of more idealized turbulent dynamo models and 
simulations (Cattaneo 1999; Cattaneo et al. 2003). In view 


of the poorly understood dependence of MHD turbulence 
on details of the problem studied, these simulations do not 
yield an unambiguous interpretation of the weak field ob¬ 
served on the Sun. With the advent and spectacular suc¬ 
cesses of realistic 3D radiative MHD simulations developed 
for the solar surface layers (Galsgaard & Nordlund 1996), 
it has become possible to study the problem numerically 
for conditions much closer to the solar case, where Vogler 
& Schiissler (2007), Pietarila Graham et al. (2010) and 
Rempel (2014) reported successful small-scale dynamo gen¬ 
eration with such simulations. Vogler & Schiissler (2007) 
and Pietarila Graham et al. (2010) investigated a range 
in magnetic Prandtl numbers and Reynolds Rm, but re¬ 
sults were still somewhat inconclusive for the combination 
of large Rm and small Pm. We here follow up on these 
results, with emphasis on the dependence on Pm at high 
resolution. 

2. Small-scale dynamos at low Pm 

A very similar dependence on the magnetic Prandtl num¬ 
ber has been observed in a variety of physical conditions 
where small-scale self-sustained magnetic field generation 
(dynamo action’) is expected. These conditions range from 
idealized models to astrophysical accretion disks and so¬ 
lar surface convection. In all these cases the dependence 
has been discovered through numerical simulations, the re¬ 
sults of which appeared to contradict basic notions of self- 
sustained field generation. 

That macroscopic behavior of MHD might depend critically 
on the magnetic Prandtl number was previously noted by 
Balbus and Hawley (1998) in the context of small-scale dy¬ 
namo action in accretion disks. An often-used model for 
small-scale dynamo action is the so-called fluctuating dy- 
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namo: a flow driven by an assumed external force acting on 
a large scale and with a random time dependence imposed 
on it. It is intended to be generically applicable to MHD tur¬ 
bulence, and believed to be sufficient to prove dynamo ac¬ 
tion (but not without controversy, see references in Iskakov 
et al. (2007)) n Numerical simulations (Schekochihin et al. 
2002) showedThowever, that the presence of self-sustained 
magnetic field generation in this model depends critically 
on the magnetic Prandtl number, with dynamo action ab¬ 
sent when Pm ^ 1 (viscosity lower than magnetic diffusiv- 
ity). The magnetic Prandtl number divides astrophysical 
systems into two very different regimes. It is usually ei¬ 
ther very small, as in the case of stellar convective cores, 
or very large, as in the interstellar or intracluster medium 
(Schekochihin et al. 2007). 

At higher Reynolds numbers or numerical resolution, 
Iskakov et al. (2007) found dynamo action below the ex¬ 
plicit numerical magnetic Prandtl number Pm = I in the 
fluctuating model, at low growth rates. The growth rate 
increases with magnetic Reynolds number, but at the low¬ 
est Pm achieved, « 0.1, it was still declining with decreas¬ 
ing Pm. The authors proposed that the results indicate a 
flattening of the growth rate to a finite positive value for 
Pm i 0. Functions that decline continuously to zero, would 
equally ht the data shown in their Fig. 1, however. Using 
a different assumption about the driving velocity held, 
Boldyrev & Cattaneo (2004) proposed that dynamo ac¬ 
tion remains possible at any low value of Pm, at sufficiently 
high Rm. Each of these models incorporates assumptions 
of uncertain physical realism about the small-scale velocity 
Held, however, and does not yield unambiguous results for 
the scaling with Pm of the minimum Rm required for dy¬ 
namo action. The existence of small-scale dynamo action 
therefore does not appear to be settled for astrophysically 
relevant magnetic Prandtl and Reynolds numbers. 

A physically realistic yet simple model that does not 
need an assumed small-scale velocity Held or forcing is the 
shearing box model for the flow in an accretion disk. Strong 
magnetorotational turbulence develops in rotating shear 
flows that are stable in the absence of a magnetic held 
(such as Keplerian shear). Lesur & Longaretti (2007) and 
Fromang et al. (2007), see also Fromang (2010, 2010b), and 
Simon (2011) studied the dependence of magnetorotational 
turbulence on Pm through simulations in which they explic¬ 
itly included viscosity and magnetic diffusivity. The results 
show that the amount of transported angular momentum 
increases with Pm- The dependence on Pm is similar to that 
found in the fluctuating dynamo simulations (Riols et al. 
2014). The auhors found that the presence or absence of a 
(weak) m ean vertical Held (a net flux crossing the disk, see 
also Sect. 5.2) is a crucial factor. 

Since most codes used in astrophysics do not include 
explicit viscosity or magnetic diffusivity, but leave this to 
the discretization errors and the stabilization algorithms 
of the code, the value of Pm that is to be associated 
with a particular simulation is not obvious a priori. The 
action of these algorithms is not equivalent to physically 
realistic diffusion coefficients. If the algorithm used for the 
stabilization term (such as a hyperdiffusivity) is of the 


^ For large-scale (global) dynamo action, the differences be¬ 
tween turbulence models and numerical results appear to be 
even more contentious. See, for example, Hughes et al. (2010) 
and references therein. 


same form for the momentum and magnetic Held vectors, 
however, one would expect its effect to correspond to a 
Prandtl number around unity. By varying the relative 
amplitude of the stabilization terms in the momentum and 
induction equations, the effect of varying the magnetic 
Prandtl number can be mimicked. 

In most simulations such adjustments are not made. 
This raises the question which value of the magnetic 
Prandtl number best characterizes the numerics. This can 
be studied by measuring effective transport coefficients 
from the simulations, for example, from the evolution of 
additional vectors for momentum and magnetic Held that 
are advected in the same way, but do not contribute to the 
equations of motion and induction (so-called ghost vectors). 
The value derived by Fromang et al. (2009) is Pmn ~ 2; sim¬ 
ilar numbers have been reported for other codes. Here Pmn 
indicates the Prandtl-like ratio of numerical diffusivities to 
keep a clear distinction from the physical magnetic Prandtl 
number Pm- This value, in the range where dynamo action 
is also observed in simulations with explicit diffusivities, ex¬ 
plains why magnetorotational turbulence has been obtained 
in most astrophysical codes. For the results reported below, 
the code includes an explicit process for modifying its effec¬ 
tive value of the numerical magnetic Prandtl number Pmn- 

The reason for the strong dependence on Pm has been 
discussed in terms of the ordering of the viscous and re¬ 
sistive length scales (cf. Moffatt (1961)). For Pm 1 , 
the viscous length scale, where the Held stretching takes 
place, is much larger than the resistive one, which then 
plays a negligible role (Batchelor 1950). The situation is 
quite different for Pm ^ 1, when the resistive scale is much 
larger than the viscous scale. As argued in Schekochihin et 
al. (2004b), in this case, small scales in the Held dissipate 
faster than they can be amplified by Held line stretching in 
the viscous flow. The authors also noted that in this inter¬ 
pretation, self-sustained dynamo action is not achievable in 
liquid metal laboratory experiments, because the required 
magnetic Reynolds numbers would be far beyond experi¬ 
mental feasibility. 

2.1. Transient behavior 

A remarkable observation that might also be relevant for 
the solar case has been made by E. Rempel et al. (2010), 
who studied the classical shearing box model of magne- 
torational dynamo action at Pm > 1 for the case where 
the mean held through the disk vanishes. Dynamo action 
saturates rapidly (a few orbital timescales) to a statisti¬ 
cally steady state, but after a finite time switches off again 
on an equally fast timescale. Since onset of dynamo ac¬ 
tion requires a finite seed Held to overcome magnetic dif¬ 
fusion, this inactive state is final. This type of behavior 
has been observed in other chaotic systems and is called 
type-II supertransient behavior. The duration of the ac¬ 
tive phase increases approximately exponentially with in¬ 
creasing Rm- The switch-off is explained as occurring when 
decay-facilitating fluctuations are by chance simultaneously 
present in all statistically independent subvolumes of the 
simulation. This occurs more readily at low spatial reso¬ 
lution or Reynolds number and explains the exponential 
dependence of the transient duration on Rm- 
Whether the supertransient phenomenon is relevant for 
stars is not clear, however. On the one hand, the steep 
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increase of the transient duration with Reynolds number 
argues against a possible relevance in environments with a 
high Reynolds number; on the other hand, some relevant 
timescales are also very long in stars. 

In the presence of a finite net flux that threads the disk, 
the supertransient effect cannot take place, since this flux 
cannot change (cf. Spruit & Uzdenksy (2005) for net flux 
in accretion disks). Provided the strength of this net flux 
exceeds the minimum required for the magnetorotational 
instability to grow, it can restart the dynamo process when 
it is on its way to decay in a supertransient accident. The 
result would be a large fluctuation instead of a switch-off. 
The key in this case is the conservation property of a net 
flux crossing a disk. The flux of an initial field consisting 
of loops that are closed within the box, for example, or a 
disk-crossing flux consisting of canceling polarities, is not 
conserved. It does not have the same effect as a net flux; it 
just acts like an arbitrary initial condition whose imprint 
vanishes with time. 


A generalization of the net flux crossing an accretion 
disk is possible. A mean field that can be considered as 
being imposed externally, for some reason, would have the 
same effect, for example, a large-scale fiel d pro duced by 
an independent dynamo process (cf. Sect. 5.21. To con¬ 
clude these considerations, the presence of some low level 
of magnetic flux from an external source may be essential 
for the long-term stability of small-scale dynamo action at 
any value of the m agne tic Prandtl number (see also the 
speculations in Sect. 5.2). 


2.2. Experimental evidence 

Liquid metals have low magnetic Prandtl numbers (on the 
order of 10“®) that are conveniently in the astrophysically 
relevant range. Reaching the high Reynolds numbers ex¬ 
pected to be needed for dynamo action has been more chal¬ 
lenging. The most successful experiment so far has been re¬ 
ported by Monchaux et al. (2007). In a turbulent shearing 
flow between counter-rotating plates in liquid sodium, these 
authors obtained dynamo action at Rm ~ 50 in the form 
of a steady field with superposed strong fluctuations. This 
success appeared to be related to a peculiarity of the exper¬ 
imental device. Dynamo action was absent until one of the 
rotating parts, made of stainless steel (low magnetic per¬ 
meability), was replaced by an iron part (high permeability 
and magnetic remanence). The remanent magnetization of 
the iron part may have played a role. Once magnetized 
by the steady component of the magnetic field, the iron 
part would have maintained a minimum field strength in 
its neighborhood. This is significant in view of the expe¬ 
rience (for example in the shearing box simulations men¬ 
tioned above) that even a weak externally imposed field 
component strongly facilitates dynamo action. Whatever 
the precise interpretation, however, the relevance of this 
experimental result for astrophysics is questionable since 
astrophysical fluids with the magnetic properties of solid 
iron are unknown. 


2.3. Situation for the Sun 

How much the dynamo mechanism depends on the details of 
the physics included is under debate ( e.g., Moll et al. (2011) 
and references therein). A range of models has been used to 


study local solar dynamo action, reaching from Bossinesq 
models (Cattaneo 1999, 2003) and compressible convective 
models in polytropic layers (Brummell et al. 2010; Bushby 
et al. 2010, 2012; Bushby & Favier 2014) to realistic so¬ 
lar surface models (Stein 2003; Vogler & Schiissler 2007; 
Pietarila Graham 2010; Rempel 2014). Bushby & Favier 
(2014) found that the critical magnetic Reynolds number 
increases with decreasing magnetic Prandtl numbers; they 
achieved a minimum value of Pmn ~ 0.1. In realistic 3D 
MHD solar surface simulations, magnetic Prandtl numbers 
Pmn ~ 1 are accessible, and for these values small-scale 
dynamo action has been found (Vogler & Schiissler 2007; 
Pietarila Graham et al. 2010; Rempel 2014). In terms of dy¬ 
namo behavior. Pm = 1 still belongs to the large magnetic 
Prandtl number limit, however. In view of the inconclusive 
results discussed above, the question is still open whether 
low- Pm small-scale dynamo action is to be expected on the 
solar surface. 


2.3.1. Indications of solar small-scale dynamo action 

Livingston & Harvey (1971) discovered an intrinsically 
weak small-scale internetwork field on the Sun. It is spread 
approximately uniformly across the solar disk and seems to 
be independent of the solar cycle. Its properties have been 
studied in detail by Martin (1988, 1990) and Zirin (1985). 
Durney et al. (1993) and Petrovay & Szakaly (1993) sug¬ 
gested that this component might be due to small-scale dy¬ 
namo action located near the solar surface. Alternatively, 
the weak-field component might represent fragments of ac¬ 
tive regions rising through the convection zone, or be a 
by-product of the decay of active regions (e.g. Spruit et al. 
(1987)). In this decay hypothesis a correlation between the 
quiet-Sun magnetic field and the solar cycle would be ex¬ 
pected. This is not evident in the observations, therefore it 
would require that the decay from the large-scale magnetic 
field to the smallest scales exceeds the solar cycle time, in 
this interpretation (Lites 2011). On the other hand, Parnell 
et al. (2009) found that the magnetic flux distribution be¬ 
tween 10^^ — 10^® Mx can be described by a single power- 
law function. This would indicate that the whole field is 
produced by the same process. Since weak fields tend to 
be compressed into strong fields by the granulation flow, 
(some fraction of) the intrinsically strong small-scale mag¬ 
netic field might be unrelated to the sunspot cycle, but 
instead result from a small-scale dynamo mechanism. The 
origin and possible variation of the strong-held component 
is of special interest because of its brightness contribution 
to the total solar irradiance (Schnerr & Spruit 2011; Foukal 
et al. 2006; Afram et al. 2011; Thaler & Spruit 2014). 


3. Calculations 

3.1. Numerical methods 

The numerical simulations were realized with the 3D 
magnetohydrodynamics code STAGGER (Nordlund & 
Galsgaard 1995; Galsgaard & Nordlund 1996). The code 
solves the time-depended magnetohydrodynamics equa¬ 
tions by a sixth-order Hnite-difference scheme using hfth- 
order interpolations for the spatial derivatives, while the 
time evolution is calculated using a third order Runge- 
Kutta scheme. For every time step the radiative trans¬ 
fer equation is solved at every grid point assuming local 
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thermal equilibrium. This is done by using a Feautrier-like 
scheme along the rays with two angles plus the verti¬ 
cal and four (j) angles horizontally, which adds up to nine 
angles in total. To incorporate the wavelength dependence 
of the absorption coefficient, the Planck function is sorted 
into four opacity bins. The equation-of-state table is cal¬ 
culated using a standard program for ionization equilibria 
and absorption coefficients (Gustafsson et al. 1973) and us¬ 
ing opacity distribution functions identical to those used 
by Gustafsson et al. (1975), which are described in more 
detail in Stein & Nordlund (1998), Nordlund (1982) and 
Nordlund & Stein (1990). For a more detailed description of 
the STAGGER code see Nordlund & Galsgaard (1995), and 
Beeck et al. (2012). The horizontal boundaries are periodic, 
while at the top and bottom are open transmitting bound¬ 
aries. The effective surface temperature resulting from the 
model is controlled through the entropy value of gas flow¬ 
ing in at the lower boundary, which is kept fixed during 
the simulation. The magnetic field is kept vertical at the 
lower boundary, thus allowing horizontal displacements of 
the field lines. At the top boundary the horizontal field com¬ 
ponents are determined from the vertical component by a 
potential held extrapolation. 

3.2. Setup 

The setup we used consists of a box with horizontal dimen¬ 
sions of 3 Mm X 3 Mm and 1.3 Mm vertical, extending 
to 475 km above the photosphere and 836 km below. The 
initial magnetic held conhguration consists of a checker¬ 
board pattern with a vertical magnetic held of alternating 
polarity and strength of 10 mG or 1 mG. This is similar to 
the setup used in Vogler & Schiissler (2007) and Pietarila 
Graham et al. (2010), who used a 4 x 4 checkerboard in a 
box with horizontal dimensions of 6 Mm x 6 Mm with a 
vertical depth of 1.4 Mm in parts of their experiments. 

3.3. Implementation of the numeric magnetic Prandtl 
parameter P^n in the STAGGER code 

In the STAGGER code the stabilization at small length 
scales is implemented through hyperdiffusivities. As already 
mentioned, the numerical magnetic Prandtl number Pin„ 
is dehned as the ratio between the magnetic hyperdiffu- 
sivity and the numerical hyperviscosity. Both of them are 
implemented following the same scheme, but by using a 
different prefactor, a variation in P^n can be produced. A 
detailed description of the diffusion scheme can be found 
in Nordlund & Galsgaard (1995) (http://www.astro.ku. 
dk/-kg/Papers/MHD_code .ps .gz) and Stein and Nordlund 
(1998). A similar procedure is used in other solar surface 
convection codes to implement hyperviscosity, which can 
be compared qualitatively but not quantitively with the 
procedure from the STAGGER code (Vogler et al. 2005; 
Pietarila Graham et al. 2010; Rempel et al. 2009; Rempel 
2014), and therefore the derived magnetic Prandtl numbers 
cannot be directly compared either. 

4. Results 

A set of simulations was made to investigate the depen¬ 
dence of magnetic held amplihcation on the numerical mag¬ 
netic Prandtl parameter Pnm- The dependence on numeri- 



Fig. 1. Evolution of (R^)^/^ as a function of the numerical 
magnetic Prandtl parameter Pmn for a numerical resolution 
of Ax = 7 km and Azq = 6 km. The solid line denotes the 
initial held strength of 10 mG, the dashed line the strength 
of 1 mG. 



X [Mm] 

Fig. 2. I Bz\ in kG at the photosphere for Pmn = 5 after 
t=231 min for a numerical resolution of Aa: = 7 km and 
Azq = 6 km. The image scale is exponential (nal°'^^) to 
amplify small structures and still be able to clearly identify 
high values 


cal resolution, as a proxy for dependence on Reynolds num¬ 
ber was investigated with an additional set of three runs. 
Figure [T] shows the dependence on Pmn- The boundary be¬ 
tween presence and absence of dynamo action appears to 
be close to Pmn = 2 . For Pmn = 5 the amplitude of the 
held increases exponentially with a growth time of ~ 200 
min, saturating at « 115 G. 

Figure shows simulation runs for an initial state of 10 
mG and a numerical Prandtl parameter Pmn = 1 for differ¬ 
ent numerical resolutions. After a short amplihcation phase 
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of about 10 min, a decline of the magnetic energy sets in. 
The code uses a uniform resolution in the two horizontal 
coordinates (Aa: in Fig. 1^; the vertical resolution is non- 
uniform. It is chosen hipest near the photosphere where 
cooling by radiation drives the flows (grid spacing at the 
photosphere is denoted by Azq in the figure). The figure 
shows that the decay time increases with increasing resolu¬ 
tion, from about 50 minutes at Aa; = 14 km to 70 min at 
7 km. At the (computationally expensive) resolution of 3.5 
km the limited time coverage of the run indicates decay at 
an even slower rate. Extrapolation to even higher resolution 
is uncertain from these data, however. 

Since the dissipation built into the code decreases with 
decreasing grid spacing, the inverse of resolution is a nu¬ 
merical analog of a Reynolds number. Since the nature of 
the numerical dissipation is rather different from physical 
dissipation and strongly depends on the order of the spa¬ 
tial discretization used, however, it is not possible to trans¬ 
late the grid spacing in our simulations meaningfully into 
effective magnetic Reynolds numbers. The uncertainty is 
smaller for the magnetic Prandtl number, since it measures 
the ratio of two quantities which, although both artificial, 
are based on the same algorithm. 



Fig. 3. Evolution of as a function of numerical 

resolution, for — 1 and an initial field strength of 10 
mG. 


5. Discussion and conclusions 

The results show the same trends as were found in the 
shearing box simulations discussed in the introduction: 
the growth rate of dynamo action increases with magnetic 
Prandtl number, and no dynamo action is detected below 
a critical value Pmn ~ 1- Eigurej^ showed that the growth 
rate appears to increase with numerical resolution, how¬ 
ever, and the corresponding critical Pmn may increase as 
well. The combination of high resolution and low Pmn re¬ 
quired to study this is numerically challenging. 

The growth rate we derived for Pmn = 5 is comparable 
with the rate reported by Pietarila Graham et al. (2010) 
for their Pmnett ~ 2. Since the numerical Prandtl parameter 
Pmn is not identical with the physical Prandtl number Pm, 
but depends on the implementation of energy dissipation 


at small scales in the code, a dependence of the effective 
Prandtl number scale on the code used was to be expected. 

Assuming a difference of a factor of a few between the 
scales, the lowest Pmnett ~ 1 reported in Pietarila Graham 
et al. (2010) would correspond to a few Pmn in our simula¬ 
tions. The fact that all results reported by Pietarila Graham 
et al. (2010) showed dynamo action is thus consistent with 
our results. 

Taken together, these results clearly show that the ques¬ 
tion of small-scale dynamo action in the surface layers of the 
Sun is still unresolved. As for the fluctuating dynamo model 
discussed in the introduction, the limit of a high magnetic 
Reynolds number (or numerical resolution) combined with 
a low Prandtl number cannot be reliably extrapolated from 
currently available results. 

5.1. Relation between intrinsically weak and strong surface 
fields 

Given that the small-scale dynamo acts and produces the 
weak-field component, might (some fraction of) the intrin¬ 
sically strong component originate from it as well? This 
seems unlikely from the results obtained with our simu¬ 
lations. Our simulation with Pmn = 5 converged to an 
average absolute vertical magnetic field strength of about 
30 G at the photosphere (after 230 min simulation time) 
and small-scale magnetic structures of surface fields of up 
to 930 G evolved (see Pig. [^, nevertheless, no magnetic 
bright points were found in bolometric intensity maps, as 
was also reported by Vogler & Schiissler (2007). The rea¬ 
son for this behavior might be that in these mixed-polarity 
simulations magnetic features of high vertical field strength 
usually are surrounded by the opposite-polarity field, and 
along with the ongoing magnetic field compression process, 
part of the magnetic flux cancels out. This means that mag¬ 
netic structures of high field strength remain too small in 
extent and short-lived to be seen as bright features. In con¬ 
trast, in unipolar simulations of the same type with the 
same average magnetic field strength, strong magnetic ele¬ 
ments were detected after about lOmin of simulation time, 
starting from a homogenous background field (Thaler & 
Spruit 2014). A small-scale surface dynamo process, if it 
takes place on the Sun, is most likely not responsible for 
most of the intrinsically strong field component. It would 
also be unable to explain small-scale field regions with a 
statistically significant dominance of one polarity over the 
other. 

5.2. Role of an imposed weak net flux 

If the experience with the shearing boxes studied in accre¬ 
tion physics is an indication, the presence of even a rela¬ 
tively weak mean flux density in the simulation may have 
a significant effect on a dynamo process. For the Sun, rem¬ 
nants of active regions spreading across the surface might 
have a similar effect of lowering the threshold for dynamo 
action. The behavior within the limit of low Pm and high 
Rm at zero mean flux would then be somewhat academic for 
the solar case. Simulations that address small-scale dynamo 
action in the presence of a dispersed strong-field component 
might then be more relevant. If it were the case that for so¬ 
lar conditions (Rm 10^, Pm ~ 10“^) dynamo behavior is 
also affected by the presence of weak fields, the strength of 
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small-scale dynamo activity might be intermittent, depend¬ 
ing on the distribution of the detritus from active regions 
rather than being present around each and every granule 
on the surface, as a local convective dynamo action would 
predict. 
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